Effect of pairing correlations on incompressibility and symmetry energy in nuclear 

matter and finite nuclei 
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The role of superfluidity in the incompressibility and in the symmetry energy is studied in nuclear 
matter and finite nuclei. Several pairing interactions are used: surface, mixed and isovector depen- 
dent. Pairing has a small effect on the nuclear matter incompressibility at saturation density, but 
the effects are significant at lower densities. The pairing effect on the centroid energy of the isoscalar 
Giant Monopole Resonance (GMR) is also evaluated for Pb and Sn isotopes by using a microscopic 
constrained-HFB approach, and found to change at most by 10 % the nucleus incompressibility Ka- 
ft is shown by using the Local Density Approximation (LDA) that most of the pairing effect on the 
GMR centroid come from the low-density nuclear surface. 

PACS numbers: 21. 10. Re, 21.65.-f, 21.60.Jz 



I. INTRODUCTION 



The nuclear incompressibility and the symmetry en- 
ergy are closely related to the isoscalar Giant Monopole 
Resonance (GMR) [l], [2| and to the isovector Giant 
Dipole Resonance (GDR) [3], respectively. The question 
of the effect of pairing correlations on the centroid energy 
of the GMR has been first addressed in Ref . [J] , and has 
recently known a renewed interest p, [6j . In both the Sn 
and Pb isotopic chains, a specific increase of the GMR 
energy, associated with the corresponding finite nucleus 
incompressibility Ka, has been predicted for the doubly 
magic 132 Sn and 208 Pb nuclei [|, 0]- A part of this ap- 
parent stiffness of doubly magic nuclei may be related to 
pairing effects which in fact decrease the GMR energy in 
open-shell nuclei. However, this study has been under- 
taken only with a pure surface pairing interaction. It is 
therefore relevant to analyze more systematically this ef- 
fect using various pairing functionals. It should be noted 
that the surface versus mixed nature of the pairing in- 
teraction is still under discussion. For instance, a recent 
systematic study based on the odd-even mass staggering 
seems to slightly favor a surface type of pairing interac- 
tion [||. 

The apparent decrease of incompressibility in super- 
fluid nuclei raises the question about a possible similar 
effect in infinite nuclear matter: until now, the nuclear 
matter incompressibility is evaluated by neglecting the 
pairing part of the functional. However, considering re- 
sults for finite nuclei, the equations of state used for neu- 
tron stars and supernovae predictions should take into ac- 
count pairing effects in the calculation of the incompress- 
ibility modulus. Therefore the question of the behavior 
of Koq with respect to the pairing gap is raised since it 
seems clear from nuclear data that the finite nucleus in- 
compressibility Ka decreases with increasing pairing gap 
[6|. A similar study for nuclear matter, as well as a more 



systematic study in finite nuclei, should be undertaken. 
This is the goal of the present work. It should be also 
noted that we will not consider the neutron-proton T=0 
pairing channel since the nuclei considered are far from 
N=Z. 

The density dependence of the symmetry energy is one 
of the most debated issues in nuclear physics at present. 
In fact, this has relevant implications (i) for nuclear struc- 
ture, since it has an important effect on the size of the 
neutron root-mean-square (r.m.s.) radius in neutron-rich 
nuclei, (ii) for nuclear reactions, e.g., in intermediate en- 
ergy heavy ion collisions where the isospin distribution 
of the reaction products is dictated by the density de- 
pendence of the symmetry energy, and obviously (iii) for 
the description of neutron stars. Review papers have 
been devoted to this topic [9|, [UJ . Empirical information 
on the symmetry energy can be obtained from various 
sources, none of them being so far conclusive by itself. 
No measurement of the neutron skin is available which is 
accurate enough to constrain the symmetry energy. The 
properties of the isovector GDR, of the low-lying elec- 
tric dipole excitations, and of the charge- exchange spin- 
dipole strength have been suggested as constraints (see 
e.g. [11(). In addition, different model analysis of heavy- 
ion collisions have been proposed as a test of the main 
trend of the symmetry energy at densities below satura- 
tion. However, in none of these studies, to our knowl- 
edge, the problem of the pairing effects on the symmetry 
energy has been addressed. 

In this work, the effects of the pairing correlations on 
the incompressibility and on the symmetry energy are 
studied consistently in nuclear matter and in finite nuclei. 
The effects coming from the correlation energy associated 
with the pairing force are included. These pairing effects 
are studied in Section II. In Section III, a Local Den- 
sity Approximation (LDA) approach to the problem is 
employed in order to understand the connection between 
the effects in infinite matter and finite systems: the 120 Sn 



nucleus is used as a benchmark. Finally, in Section IV, a 
microscopic study on the role of superfluidity in the in- 
compressibility of finite nuclei is undertaken, employing 
several pairing interactions: surface, mixed and isovector 
dependent. 



II. NUCLEAR MATTER 

In this Section, we study the effects of the pairing cor- 
relations on the incompressibility and the symmetry en- 
ergy in nuclear matter. 

A. Energy density 



TABLE I: Strength v (in MeV-fm 3 ) of the pairing interac- 
tions obtained in the case of various Skyrme functionals. The 
values of vo are adjusted, for each 77, to obtain equivalent 
results to those of Ref. [Q for the two neutron separation 
energy in the Sn isotopes by the HFB calculations with the 
parameter set SLy5. The energy cutoff for the pairing window 
is taken to be 60 MeV. 





77=0.35 


77=0.65 


77=1.00 


SLy5 [12] 


-285 


-390 


-670 


Sk255 [14] 


-265 


-390 


-600 


Sk272 [14] 


-265 


-390 


-600 


LNS [15] 


-250 


-390 


-670 



The nuclear energy density e (e/p = E/V) is the sum 
of the Skyrme part, esicyrme, that includes the kinetic 
energy [12j, plus the pairing energy density, 



Here 



eSkyr 



e pair = --(N n A 2 n + N p A 2 p ) 
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In Eq. (J2), A T is the pairing gap and N T is the density of 
states, given by N T = m*fci? T /7r 2 /i 2 , with r = n,p. The 
energy density e is a function of the total density p and of 
the asymmetry parameter S = (p n — p P )/p- In the T = 1 
channel, several pairing interactions are defined by 
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as a function of the value of 77 that can range from 
(volume- type pairing) to 1 (surface- type pairing). In Eq. 
([3]) the parameter a is set to 1 and po is taken as the sat- 
uration density of symmetric nuclear matter throughout 
all the study; moreover, we adopt the parameters 77=0.35 
and 0.65 for the volume-surface mixed-type pairing inter- 
actions, and ?7=1.0 for the surface-type interaction. The 
values of vq in all these cases are adjusted, for each 77, 
in such a way to obtain equivalent results for the two 
neutron separation energy in the Sn isotopes by HFB 
calculations with the SLy5 parameter set. The pairing 
cutoff energy is set at 60 MeV [1J|. These values of tj 
are given in Table HI In the following, these pairing in- 
teractions will be denoted as IS, because they depend on 
the isoscalar density. 

We have also considered pairing interactions having 
the isovector density dependence with 8 = (p n — p p )/p in 
addition to the isoscalar density dependence. The MSH 
interaction is defined as Hal 
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with v = -448 MeVfm 3 , 77 s = 0.598, a s = 0.551, r) n = 
0.947, a„ = 0.554 (with a cutoff energy of 60 MeV). The 
YS interaction has also the isospin dependence as [13] 
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with v Q = -344 MeVfm 3 , 770 = 0.5, 771 = 0.2, 772 = 
2.5 (with a cutoff energy of 50 MeV). The parameters 
mentioned have been used in connection with the SLy5 
interaction. In the following these pairing interactions 
will be denoted as IS+IV. 

Effective interactions in the pairing channel are faced 
to double counting problem [18|. However using a zero 
range interaction constrained on bare interactions allows 
to avoid this problem (l9j . This motivates the use of a 
different interaction in the the pairing channel compared 
to the particle-hole one. The pairing interactions used in 
this work are of two types: either fitted on BCS gaps in 
symmetric and neutron matter calculated with bare in- 
teraction, as in the MSH case, or designed to fit observ- 
ables in nuclei, as in the IS and YS cases. These pairing 
interactions are featured with zero range and with a cut- 
off, following the prescription of Ref. 19]: our aim is to 
provide reasonable pairing description in nuclei in order 
to use the same interaction in nuclei and in nuclear mat- 
ter. The MSH interaction is considered as an extension 
of Ref. [19| including isospin dependence. It has been 
shown that it is possible to study the surface properties 
of the pairing interaction using slabs of nuclear matter 
[20T [21| . However, it should be mentioned that one usu- 
ally deals with only the first order in the diagrammatic 
expansion of the many body equations. Thus adjusting 
the pairing gap on the bare interaction is done as this 
level and is therefore perfectible: there are screening ef- 
fects which are at next order for instance. Therefore ad- 
justing the gaps on the bare interaction is complementary 
to other ways such as constraining the pairing interaction 
to fit the gaps in finite nuclei. We use both approaches 
as explained above. 



In the present work, symmetric nuclear matter is stud- 
ied as well as the behavior of the symmetry energy. On 
this purpose, it should be noted that the MSH interac- 
tion is obtained by constraining the neutron and proton 
gaps of the bare interaction in both symmetric and pure 
neutron matter. The YS interaction is adjusted on sev- 
eral nuclei with different isospin and there is an explicit 
isospin dependence. To perform a study in asymmetric 
nuclear matter, the non linearity of the energy gap on 
the isospin degree of freedom should be considered [22[ 
and may be studied in a forthcoming study. 

The pairing gap in uniform matter is obtained from 
the BCS gap equation [23| 



A* 



£ 



-Vkk< 



2E k , 



(6) 



solved under the condition of the particle number con- 
servation. In a given volume V one assumes constant 
density given by 
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where the quasiparticle energy is defined as Ek >T 

\/{ e r{k) 



A 2 fc , r 



p T )^ -+- la 1 ^ , e T (k) is the single particle en- 
ergy, and /i T is the chemical potential. In Eq. (J5J), Vkk' is 
the pairing matrix element for the plane waves, namely 
(kk\v\k'k') . Notice that in the case of the zero-range 
pairing interaction, the pairing gap A& is independent of 
k. 

In Fig. [T] we display the pairing gap A T , the pairing en- 
ergy per particle e pa ir/p and the percentage of the pairing 
energy with respect to the total energy e in symmetric 
matter for the various pairing interactions together with 
the SLy5 Skyrme interaction in the mean field channel. 
There is a critical density p c sw 0.11 fm -3 at which all the 
pairing interactions give almost the same result for the 
pairing gap. This has already been noticed in Ref. [131 ] 
and may be related to the fact that in fitting the two- 
neutron separation energy one is sensitive to the space 
region of the nuclear surface, where the density is some- 
what lower than the saturation density: therefore the 
pairing gap is constrained rather at p c than at p$. No- 
tice that the parameters of the MSH pairing interaction 
have been determined without using constraints from fi- 
nite systems, and A at po does not necessarily coincides 
with the one from the other pairing interactions. Above 
p c , the more surface- type the pairing interaction (that 
is, the larger r\ is taken), the smaller the pairing gap 
A r . Below the critical density, the trend is reversed: the 
more surface-type the pairing interaction is, the larger 
the pairing gap. The contribution of the pairing energy 
is increased at low densities. Around the saturation den- 
sity, the pairing energy per particle is much smaller than 
the binding energy (—16 MeV). 



B. Incompressibility and symmetry energy 



The compressibility \ is usually defined by 
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and the pressure is related to the energy density e by 
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The chemical potential is defined by 
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,/From Eqs. ([5]) and ©, we obtain 
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The compressibility x{p) is a function of the density p and 
the asymmetry parameter 5. Furthermore one defines the 
incompressibility at the saturation density in symmetric 
nuclear matter by 
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The relation between K^ and x is given by 
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It is worth keeping in mind that the incompressibility K x 
is defined only at the saturation density po, and in par- 
ticular the relation (|13[) is valid only for p = po- Indeed, 
more generally, we can define the density-dependent in- 
compressibility (or bulk modulus) as [24l j 
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which coincides with the incompressibility K^ at the sat- 
uration density. 

The symmetry energy S(p) is defined by 
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1 d 2 e/p 



2 dd 2 
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The symmetry energy can be expanded, around the sat- 
uration density, as 



S(p) = J + L 



P~ Po 
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where J is defined by J = S(po), L — 3po^f \p , and 
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FIG. 1: (Color online) Pairing gap A T , pairing energy per particle e pa ir/p, percentage of pairing energy with respect to the 
total energy e and equation of state around the saturation point in symmetric matter obtained with various pairing interactions 
employed in connection with the SLy5 Skyrme interaction. The solid, dotted and dashed lines correspond to the IS pairing 
interactions with r\ = 0.35,0.65 and 1.0, respectively, in Eq. (3). The dashed-dotted and dashed-dotted-dotted lines show the 
results of IS+IV interactions in Ref. [15] (MSH) and Ref. [16] (YS), respectively. See the text for details. 



C. Results 

Fig. [5] displays the pressure P in Eq. ©, the incom- 
pressibility K{p) in Eq. (fT4"|) . and the symmetry energy 
S(p) in Eq. (TTB")) without pairing (top panels), and the 
contribution of pairing to these quantities (bottom pan- 
els), using the SLy5 interaction. This contribution is 
calculated with the same equations, but considering only 
the pairing term of the energy density in Fig. [TJ The 
same pairing interactions have been considered here as 
in Fig. [TJ Close to the saturation density, the contribu- 
tion from pairing is very small. This is also illustrated 
in Table [TTJ the pairing interaction has small effects at 
the saturation density. In the case of the incompress- 
ibility K^, pairing can still produce a few % effect (for 
instance K^, is changed from 230.2 MeV to 223.9 MeV in 
the case of the MSH pairing interaction). The MSH, YS, 
IS 0.35 pairing interactions modify the incompressibility 
by 3 to 6 MeV, that is, by about 2%. It should be noted 
that at the saturation density, the contribution to the 
slope parameters of the symmetry energy, L, and K sym , 
of the interactions MSH and YS is larger than that of 
the other IS forces. The effects on L can be about 15% 
while K sym can be modified in an important way. This 
is related to the dependence of these pairing interaction 
on the isovector density. 

However at lower densities, the pairing effects become 
appreciably larger as seen in Fig. [2j In the case of 
the pure surface pairing, there are important contribu- 
tions to the pressure, incompressibility and symmetry 
energy: these quantities can be strongly affected by pair- 
ing, which can lead to variations up to about a factor 2. 
Other pairing interactions also provide significant correc- 
tions to the pressure and the incompressibility, typically, 
around 10%. In the case of the symmetry energy, below 



TABLE II: Properties of nuclear matter for various pairing 
interactions. The SLy5 Skyrme force is used for the mean 
field. 



Pairing po 
[fm" 3 ] 


E/A(p ) 

[MeV] 


[MeV] 


J 

[MeV] 


L 

[MeV] 


-"-sym 

[MeV] 


no pairing 0.1604 


-15.999 


230.2 


32.03 


48.25 


-112.3 



IS 77=0.35 0.1601 -15.998 227.3 31.93 48.49 -129.7 

IS 77=0.65 0.1603 -15.998 228.1 32.02 48.30 -113.7 

IS 77=1.00 0.1604 -15.999 230.1 32.03 48.25 -112.3 

MSH 0.1599 -15.998 223.9 31.33 55.77 -139.7 

YS 0.1602 -15.998 227.0 31.39 52.04 13.2 



p Ri 0.1 fin -3 , the IS+IV pairing interactions (MSH and 
YS) predict an opposite and positive contribution com- 
pared to the negative contributions of IS pairing interac- 
tions. It should be noted that the pairing contribution 
to these quantities is generally larger at densities below 
saturation. 

To obtain a more general view of the pairing effect on 
the incompressibility, Table QTI] displays the K^ values 
obtained for SLy5, LNS, Sk255 and Sk272 Skyrme func- 
tionals, with various pairing interactions. The correction 
induced by the pairing interaction IS 0.35 is largest one 
among the IS interactions and reduces the incompressibil- 
ity Koo by about 3 MeV. The MSH interaction induces 
a correction of 6.3 MeV on the incompressibility (Table 
HI)) . It should be noted that the pure surface pairing in- 
teraction provides no modification of Koo. Depending on 
the Skyrme models, there shall also be an effect due to 
the different effective masses m* /m, but they are incor- 
porated in the renormalization of the pairing interaction 
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FIG. 2: (Color online) Pressure, incompressibility, and symmetry energy without pairing (top panels), and the contributions 
of various pairing interactions to these quantities (bottom panels). The HF energy eskyrme is calculated using the SLy5 
interaction. For details, see the caption to Fig. 1 and the text. 



TABLE III: Nuclear matter incompressibility K x (MeV) for 
SLy5 QJ], LNS 0, Sk255 QJ] and Sk272 [3 Skyrme func- 
tional. The dependence of Koo on the pairing interaction is 
displayed: mixed (IS r7=0.35), surface (IS 77=1. 00). 



Pairing 


SLy5 


LNS 


Sk255 


Sk272 


no pairing 


230.2 


211.0 


255.2 


271.8 


IS ?7=0.35 

IS 77=1.00 


227.3 
230.1 


208.4 
211.0 


251.3 
255.2 


268.3 
271.8 



parameter vq shown in Table [I] 

It is expected that the above pairing effects at low den- 
sities may also affect finite nuclei. In the case of incom- 
pressibility, we can define a finite nucleus value Ka and 
expect that this value is affected by the pairing more than 
Koo, due to the presence of a lower density region, i.e. 
the nuclear surface. We analyze this point in the next 
Section, and we argue that a similar reasoning holds for 
the symmetry energy. 



III. LOCAL DENSITY APPROXIMATION 

This section relates the general expressions in uniform 
matter obtained in Sec. [IT] with the observables in finite 
nuclei in the local density approximation. The aim is to 



estimate the role of pairing in the incompressibility and 
symmetry energy of finite nuclei in a simple and trans- 
parent way. The validity of the LDA will be estimated 
by comparing the predicted nuclei incompressibility with 
the one obtained by a microscopic approach. 

The binding energy per nucleon in the LDA reads 



B A (N,Z) = B mcl (N,Z)+B Co 



'A 4 /3 



(17) 



where -Bnuci(-^, Z) includes the bulk, surface and pairing 
contributions. It is defined by 



B Nucl (N,Z) 



d A r e(r) 



(18) 



where e(r) = e(p n (r),p p (r)) = e Skyrme (r) + e pair (r) as 
was defined in Eq. (1). The neutron and proton densities 
p n (r),pp(r) can be obtained, in the present context, by 
means of a spherical HF calculation. The pairing contri- 
bution to the binding energy is defined by 



B pail .(N,Z) = jJd 3 



c pair 



(r). 



(19) 



-BnucI can be expanded around the saturation density, 
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FIG. 3: (Color online) Radial dependence of e(r), Kj^ uc i,(r) and SU(r) for the various pairing interactions considered (bottom 
panel) using the SLy5 force in 120 Sn within the LDA. The contribution of the Skyrme term plus the kinetic term is displayed 
in the top panel. See the captions to Fig. 1 and the text for details. 



where the incompressibility in nuclei, Ka, is defined by 
K A = i^Nud + i^coui • Z 2 A- A /\ and 



2 d 2 B Nucl (N,Z) 

-ftNucl = 9Po a~2 : 



(21) 



while the pairing contribution to the incompressibility is 
defined by 



-t^pair — ^^1 



2 d 2 B pair (N,Z) 
dp 2 



(22) 



The Coulomb contribution, i^coui, can be evaluated us- 
ing for instance the Thomas-Fermi approximation [cf. 
Eq. (Al) in Ref. 25]]. It will not be included in the 
present work but the value obtained in Ref. 25] is —8 
MeV < Kcou\ < — 4 MeV, depending on the interaction 
which is used. 

The symmetry energy in finite nuclei, Sa, is defined 
by 



S A = 



ld 2 B Nucl (N,Z) 



dS 2 



8=0 



and the contribution of the pairing correlations to Sa is 
defined by 



_ l d 2 B pair (N,Z) 

■Vir - 2 gp 



(24) 



Introducing the mass formula (fT8"f into Eq. (|2T|) . one 
obtains 
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Nucl 



.4 



d 3 r K Nuc i(r) 



with 
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Pa 
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For small values of the density (p < 0-6po, that is r 
> 5 fm in 120 Sn), the incompressibility is found to be 
negative: this is due to the spinodal instability in nuclear 
matter which is not present in finite systems [26] . For this 
reason, the integral (|25p is limited to the region where 
-Knuci(?~) is positive. In this way, the spurious component 
due to the spinodal instability is removed. 

Introducing the quantity 



o , s 1 d 2 e 
SA{r) = 2^W 



= -S( P (r)), 

8=0 P0 



(23) the symmetry energy in nuclei (|2"3")l reads 



c P0 



d 3 r S A (r). 
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5=0 



We first perform a self-consistent HF calculation which 
provides the neutron and proton densities in 120 Sn. From 
these densities we deduce the radial distributions of mean 



TABLE IV: Contributions of pairing correlations to the bind- 
ing energy, the incompressibility and the symmetry energy in 
120 Sn. The mean field is calculated by using SLy5 interaction 



Pairing 


-Dpair 


-"-pair 


•Jpair 




[MeV] 


[MeV] 


[MeV] 


IS 77=0.35 


-0.03 


-0.5 


-0.25 


IS 77=0.65 


-0.11 


-3.9 


-0.03 


is 77=1.00 


-0.64 


-6.0 


-0.03 


MSH 


-0.13 


-3.2 


-0.93 


YS 


-0.05 


-1.9 


-0.24 



field part and pairing part of e(r) given in Eq. (1) , 
K^ uc i(r) given in Eq. ([2"6"j) . and Sa(t) given in Eq. (1271) : 
these radial functions are shown in Fig. [3] As expected 
from the results discussed in the previous section, the 
pairing effects on e(r), K^ uc \(r) and Sa{t) come from 
the low density surface region. 

From Eqs. ([15]). (j2"5|) and (J25J), we obtain, in the 
SLy5 case, B A = -13.5 MeV, if Nud =119.8 MeV, and 
SU=25.7 MeV without the contribution due to the pair- 
ing correlations. The Coulomb contribution has not 
been included. The value for -KnucI should be compared 
with that of 141 MeV obtained by the constrained HFB 
(CHFB) calculations presented in the next section (cf. 
Fig. |BJ). It should be noted that the CHFB calcula- 
tions take into account the contribution coming from the 
Coulomb interaction. This contribution is estimated to 
be about 20 MeV in 12 "Sn. using the values oiKcoul from 
Ref. [25J . The good agreement between the LDA and the 
CHFB results ensures that LDA provides a sound frame- 
work to relate the nuclear matter incompressibility and 
the finite nucleus one. 

The contributions of pairing correlations to the bind- 
ing energy, the bulk modulus and the symmetry energy 
are shown in Tab. IIVI for the various pairing interactions 
considered. The contribution of the surface-type pairing 
(IS 77=1.0) reduces Ka by about 5%, whereas, for the IS 
mixed-type (r7=0.35 or 0.65) and the IS+IV (MSH and 
YS) pairing interactions, the effect on Ka is predicted to 
be smaller. In Table IIV| it is also observed that pairing 
effects affect the binding energy by few percents, up to 
5% for the surface-type pairing interaction. For the sym- 
metry energy, pairing effects are negligible, being below 
1% except the IS+IV pairing (MSH). 



IV. FINITE NUCLEI 

In the previous section, the LDA has shown that 
the pairing effect on the symmetry energy is negligible, 
whereas appreciable effects are observed in the case of 
the incompressibility. In this section, the role of pair- 
ing effects on the finite nucleus incompressibility Ka is 
investigated using a microscopic approach. 



We use the sum rule approach in order to calculate the 
centroid energy of the isoscalar GMR. It is known that 
the finite nucleus incompressibility Ka is related to that 
centroid energy by means of the relation 



£1 



ISGMR 



IH 2 Ka 

m(r 2 ) 



(29) 



where m is the nucleon mass and (r 2 ) denotes the ground- 
state expectation value of the square radius. In a micro- 
scopic approach, for the so-called scaling Ka, we calcu- 
late the energy as 



-Eisgmr = 



mi 

TO_l 



(30) 



where the k — th energy weighted sum rule is defined as 

m fe = 5>*|(*|Q|0}| 2 , (31) 

i 

with the RPA excitation energy Ei and the isoscalar 
monopole transition operator, 



E 



r?. 



(32) 



The m\ moment is evaluated by the double commutator 
using the Thouless theorem [271 ]: 



mi 



2ft2 / 2\ 



(33) 



In the present HFB calculations, the energy cutoff is 60 
MeV, and the j max value is 15/2 in the case of IS pairing 
and extended to 25/2 for the IS+IV MSH pairing, which 
are the cutoffs used in the design of these pairing inter- 
actions [1J, [l6|], in order to ensure convergence of the 
results. It should be noted that the strength v of the 
IS pairing interactions is adjusted in nuclei for its corre- 
sponding jmax and energy cutoff. Therefore the different 
values of j max between IS and IS+IV pairing has little 
influence on the calculations. 

Concerning the evaluation of the m_i moment, the 
CHFB approach is used. It should be noted that the 
extension of the constrained HF method [2|, [28| to the 
CHFB case has been recently demonstrated in Ref. [29|] 
and employed also in Q . The CHFB Hamiltonian is built 
by adding the constraint associated with the IS monopole 
operator, namely 



H,, 



= H + \Q, 



(34) 



and the m_i moment is obtained from the derivative of 
the expectation value of the monopole operator on the 
CHFB solution |A), 



1 

m_i = 

2 



|<AMA) 



(35) 
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FIG. 4: K x versus K A for 208 Pb obtained by the CHF 
method for several Skyrme interactions. 
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FIG. 5: K x versus K A for 114 Sn obtained by the CHF and 
the CHFB method with surface-type and mixed-type pairing 
interactions for several Skyrme interactions. 



We first investigate the correlations between nuclear 
matter incompressibility K^ and finite nucleus, Ka, in 
the absence of pairing correlations. This correlation has 
been found in all the previous papers on the subject, but 
is reported here as a benchmark for further considera- 
tions concerning the effect of pairing. For this purpose 
we choose the doubly magic 208 Pb nucleus. Figure |4] dis- 
plays Ka (obtained with the CHF method) versus K x 
for the four Skyrme interactions LNS, SLy5, Sk255 and 
Sk272. These four interactions span a large range of in- 
compressibilities and have been fitted by using different 
physics inputs: the neutron matter EOS from realistic 
forces in the case of SLy5, Briickner-HF calculations in 
nuclear matter in the case of LNS, and the empirical 
properties of symmetric uniform matter plus a few bind- 
ing energies and charge radii of selected nuclei (the same 
which had been used to fit some relativistic functionals 
like NL3) in the case of Sk255 and Sk272. In this sense, 
these interactions provide representative samples of the 
Skyrme functionals. 

It should be stressed that both Ka and K^ are here 
evaluated consistently in a microscopic model. As ex- 
pected, Ka is clearly correlated with K^. To find out 
the pairing effects, we also study the case of the open- 
shell nuclei 114 Sn and 120 Sn in Figs. [S]and[ni respectively. 
For each Skyrme interaction, three results are shown: (i) 
the CHF result, without pairing, (ii) the CHFB result 
using the surface-type pairing interaction, and (iii) the 
CHFB result using the mixed-type (?7=0.35) pairing in- 
teraction. In the case of SLy5, the IS+IV MSH pairing 
interaction is also used. The surface-type interaction de- 
creases the finite nucleus incompressibility Ka by about 
10% in 114 Sn and 5% in 120 Sn whereas the mixed-type 
pairing interaction has a negligible effect on Ka- Con- 
versely, it should be reminded that the mixed pairing 
interaction has some effect on K^ whereas the pure sur- 
face pairing interaction has a negligible effect on K^, 
as seen on Fig. \5\ In the case of LNS, the reduction 



of Ka is smaller for the surface pairing, and an increase 
of Ka is even observed for the mixed pairing case. In 
the SLy5 case, predictions using the IS+IV MSH pairing 
interactions show no variation of Ka but affects K x . 

To study how these conclusions in 114 Sn are sensitive 
to the nuclear shell structure, the results for 120 Sn are 
displayed in Fig. [51 where pairing effects are expected 
to be smaller than in 114 Sn due to the subshell closure. 
In this case the reduction of Ka due to the surface-type 
pairing effect drops to 5%. For the mixed pairing inter- 
action a small increase of Ka is observed. This feature 
is again more pronounced in the LNS case. In the SLy5 
case, predictions using the IS+IV MSH pairing interac- 
tions show no variation of Ka but affects K^ . It should 
be noted that similar trends are observed with the LDA 
predictions (See Table HV|). Also, they are consistent with 
previous studies |a,|g]. To further study shell effects, the 
same calculations have been performed on 126 Sn, and a 
similar pattern than Fig. [5] is found, showing that the 
present results a rather independent from shell effects in 
open shell nuclei. 

Evidently, the pairing effects tend to decrease the fi- 
nite nucleus incompressibility Ka in the surface pairing 
case, whereas K^ is decreased in the mixed pairing case. 
Hence the question of constraining the pairing interac- 
tion through precise (typical resolution of few hundreds 
of keV) GMR measurement is raised. The systematic 
softness of open-shell Sn and Cd isotopes measured re- 
cently [3CH32I I through the energy of the GMR might be 
the sign of a surface pairing interaction at work [6j . 



V. CONCLUSIONS 

The effect of superfluidity on the incompressibility has 
been studied in both nuclear matter and finite nuclei, us- 
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FIG. 6: Same as Fig. 5, but for 120 Sn. 



ing various pairing functionals. A small effect is observed 
on the nuclear matter incompressibility at the saturation 
density and the symmetry energy, but non-negligible on 



sym- 



L and K ; 

However at lower density, the pairing effect on the in- 
compressibility is significant and can have an substantial 
impact on neutron stars studies or on the interpretation 
of multifragmentation data. It has been shown that the 
LDA provides a relevant framework for a qualitative un- 
derstanding and interpretation of the microscopic results. 
The effect of the pairing correlations is localized near the 



surface of nuclei and the effect of the pairing correlations 
is to make slightly softer nuclear EOS. Especially in the 
low density region in nuclear matter, the pairing effect 
is more noticeable. This may explain why such effects 
are expected to happen in the surface of the finite nuclei. 
In the case of the IS+IV pairing interaction, no strong 
effect is observed on Ka- In general, the pairing effects 
on the finite nucleus incompressibility Ka are more im- 
portant when the interaction is more surface type (larger 
t] value). 

This study shows that with respect to current exper- 
imental uncertainties, the pairing effects should be con- 
sidered when extracting the incompressibility value from 
GMR data which can now reach an accuracy of several 
hundreds of keV [30| • Experimentally it would be useful 
to measure the GMR on isotopic chains, including both 
open-shell and doubly magic nuclei such as 132 Sn. Such 
measurements are starting to be undertaken [30l432| and 
will be extended to unstable nuclei 13311. 
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